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Abstract 

We study the behavior of a quantum particle confined to a hard-wall strip of 
a constant width in which there is a finite number N of point perturbations. 
Constructing the resolvent of the corresponding Hamiltonian by means of 
Krein's formula, we analyze its spectral and scattering properties. The bound 
state-problem is analogous to that of point interactions in the plane: since 
a two-dimensional point interaction is never repulsive, there are m discrete 
eigenvalues, 1 < m < N , the lowest of which is nondegenerate. On the 
other hand, due to the presence of the boundary the point interactions give 
rise to infinite series of resonances; if the coupling is weak they approach the 
thresholds of higher transverse modes. We derive also spectral and scattering 
properties for point perturbations in several related models: a cylindrical 
surface, both of a finite and infinite heigth, threaded by a magnetic flux, 
and a straight strip which supports a potential independent of the transverse 
coordinate. As for strips with an infinite number of point perturbations, we 
restrict ourselves to the situation when the latter are arranged periodically; 
we show that in distinction to the case of a point -perturbation array in the 
plane, the spectrum may exhibit any finite number of gaps. Finally, we study 
numerically conductance fluctuations in case of random point perturbations. 



Contents 



1 Introductiorj 



2 A one— center perturbation in an infinite strip 



2.1 Boundary conditions] 4 

272 The resolvent) 5 



2.3 Scaling behavioi] 6 



2.4 The resolvent a la Zorbas| 7 



1 



2.5 Spectral properties! 7 

2.6 Scattering! 10 

2~7 Resonances! 12 



A finite number of point interactions! 14 

O The resolvent! 14 



3.2 The discrete spectrum! 15 



3.3 Embedded eigenvalues] 17 



3.4 The limits of strong and weak coupling] 19 



3.5 More about bound states 21 



3~6 Scattering! 23 



4 A cylindrical strip in axial magnetic field 



24 

41 The free Hamiltonian! 25 



4.2 The perturbed resolvent! 26 



4.3 Spectral properties! 26 



Point perturbations on an infinite cylinder 

5.1 The resolvent and scattering amplitudes . 



27 

27 



6 An infinite strip in external field 



28 

6.1 The unperturbed Hamiltonian] 28 



6.2 Point perturbations] 29 



6.3 Scattering] 30 



7 Periodic perturbations 






7.1 The Floquet-Bloch analysis . 


7.2 A single array of perturbations 





31 



5 Random point perturbations 



36 



1 Introduction 

Point interactions were introduced in early days of quantum mechanics by R. Kronig, 
W.G. Penney, E. Fermi and others. It lasted several decades, however, before it 
became clear how they can be treated in a mathematically sound way and in what 
sense they correspond to the intuitive notion of a 5-shaped potential. This in turn 
made it possible to develop this method into a powerful tool for analyzing quantum 
systems in which the interaction is supported by a family of sets which are small on 
an appropriate scale and well separated. 

A thorough and extensive discussion of the point interaction method can be 
found in the monograph 0; references to some recent work are given, e.g., in P, |§. 
This does not exhaust, however, the list of possible applications. In this paper we are 
going to study another system with point interactions, namely a quantum particle 
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living in a strip of a fixed width d with the Dirichlet (i.e., hard wall) boundary, in 
which there is a finite number of pointlike perturbations. 

This problem has an obvious physical motivation. The strip is a natural model 
for quantum wires which are an object of intensive experimental and theoretical 
interest these days. Though the presently available technologies make it possible to 
fabricate mesoscopic structures of extremely pure semiconductor material in which 
the electron motion is predominantly ballistic, a real wire is always non-ideal, con- 
taining impurities in the crystallic lattice which affect its conductivity properties. 

Even if one treats the one-electron problem within a suitable range of energy in 
the approximation of a particle with the effective mass whose motion is confined to 
a strip, the scattering problem in the presence of impurities is not easy to solve. The 
assumption that the impurities have a point character can turn the solution into a 
more manageable task; at the same time it is a reasonable approximation because a 
typical impurity consists of an alien atom in the lattice while the wire cross section 
includes at least several but mostly many atomic layers. Moreover, an experimental 
way to produce artificial impurities has been reported plf ; though the latter are 



not exactly pointlike, their existence stresses the need for better understanding to 
scattering in quantum waveguides. 

Importance of point interactions as models of impurity scattering in a two- 



dimensional electron gas was realized long time ago EI, however, it seems to at- 
tract a new attention recently when several studies of point perturbations on strips, 
cylinders or tori has appeared — see, e.g., |^, [£| and references therein. In par- 
ticular, P. Bagwell has pointed out the role played by quasibound states at excited 
transverse modes which are unstable due to intermode coupling. Let us remark that 
similar resonances appear generically if there is a locally attractive interaction in 
the wire (the latter may be of a various origin; it can be induced also in a purely 
geometrical way — see, e.g., |TU|, [IB], [II], [H]). The most remarkable manifestation 
of the resonances are local deformations of the ideal steplike conductance pattern. 

On the other hand, the point scatterers in the above mentioned papers are 
treated in a simplified way, usually with the model space restricted to a finite number 
of transverse modes. This does not allow us to introduce the interaction in a proper 
local way. In fact, the authors of Ref.|| discuss the difference between the situation 
when the scatterers are attractive or repulsive. It is known, however, that — in 
distinction to one and three dimensions - - a two-dimensional point interaction 
obtained as a limit of a family of scaled potentials is never repulsive; this fact 
follows from the logarithmic character of singularity of the corresponding resolvent 
kernel — cf. @, || or |, Sec.1.5]. 



This gives a motivation to analyze the problem of point perturbations in a strip 
more carefully; we follow the ideas which have yielded the one-center scattering ma- 
trix discussed in [[TJJ . As in the case without the boundaries |I| , the point perturba- 
tions are constructed as self-adjoint extensions of a symmetric operator obtained by 
restriction of the free Hamiltonian in the strip. Our main tool is the Krein's formula 
which allows us to write down the resolvent of the Hamiltonian; the difference from 
the free resolvent is a rank N operator, where iV is the number of impurities. Since 
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we are interested in local perturbations only, the number of parameters is again N ; 
they are naturally interpreted as the corresponding coupling constants. 

The spectral problem is then reduced to finding eigenvalues of an NxN matrix. 
In distinction to the case of point interactions in a plane, however, the matrix 
elements are not elementary functions but sums of certain series. Nevertheless, 
one can say a lot about the spectrum. In particular, the mentioned fact about 
attractivity of two-dimensional point perturbations is not affected by the presence 
of the strip boundary. One can also derive the S-matrix which yields conductance 
modulations due to the impurities, as well as infinite series of resonances associated 
with each point perturbation. Another feature observed in scattering are turbulent 
patterns of the probability flow caused by the impurities, which have been noticed 
in Ref.|J and in || |31| in a bent-guide context. 

Furthermore, strips with point perturbations offer a useful test field for some 
properties dependent on the dimension of the configuration space, such as number 
of open gaps in a periodic situation, or existence of localization in case of random 
perturbations, because, roughly speaking, they are "halfway" between the one- 
dimensional and two-dimensional situation. We shall show that while a straight 
polymer in the plane has at most one gap |T], Sec. III. 4], a "coated polymer", i.e., a 
periodic array of point perturbations in a strip can exhibit any finite number of open 
gaps. This raises the question, of course, whether the Bethe-Sommerfeld conjecture 
can be proven in this situation. 

Finally, if the perturbation positions or coupling constants are distributed at 
random, we can study statistical properties of observable quantities such as conduc- 
tance of a wire with point perturbations. The explicit solvability of our model makes 
it suitable for comparison with predictions of random-matrix theory; this will be 
done in the concluding section. 

2 A one— center perturbation in an infinite strip 

Consider a straight planar strip Q := M x [0, d] and suppose that the free motion in 
Q is governed by the Dirichlet Laplacian [^5|, Sec. XIII. 15]; since the boundary of Q 
is smooth, it acts as H ip = —Aip on the domain of all if) for which the rhs (in the 
sense of distributions) belongs to L 2 {VL) and which satisfy the boundary conditions 



for all x G 1R . Unless stated otherwise, we put for simplicity d = ir ; the results for 
a general strip width are then obtained by a simple scaling transformation. 

2.1 Boundary conditions 

Our first goal is to construct a one-center singular perturbation to H supported by 
a point a := (a, b) with b G (0, 7r) . This can be done using the procedure mentioned 
in the introduction; since locally there is no difference between point interactions 



i/;(x,0) = i/;(x,d) = 
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in the plane and in the strip, we can use the standard boundary conditions |], 
Chap. 1. 5]. We introduce the regularized boundary values 

ip (x) 

L (ip,a) := lim - — — — — - , Li(ip,a) := lim [ip(x) — Lq(iP, a) In \x— a\] 
|af— fi\— >o \ia\x — a\ \x—s\^o 

and require 

LiO,a) + 27raL O,a) = 0. (2.2) 
To any a G M we define the self-adjoint operator H(a, a) as 

H(a,a)ip '■= — , 

(2.3) 

D(H(a,a)) := : -A-0 is L 2 on fi\{a} and (2.1), (U) are satisfied } ; 

here again — Aip is understood in the sense of distributions. The case a = oo , 
i.e. L (ifj, a) — , corresponds to the free Hamiltonian H . 

2.2 The resolvent 

To find the resolvent of H(a, a) we start with that of H . We use the decomposition 

into transverse modes, L 2 (Q) = 0^L X L 2 (M) (g> {Xn} , where Xn(y) '■— sin(ny) ; 
then the free Hamiltonian can be written as 

oo £. 

H = h n ® I n , /i n := - — + (2.4) 

n=l dX 

with D(h n ) := AC 2 (M) . It follows that the free resolvent is an integral operator 
with the kernel 

^ oo e ik n (z)\xi— xa| 

Go(xi,f 2 ; 2) = (-Ho- ^) _1 (fi,x 2 ) = - V p— sin(nyi) sin(m/ 2 ) , (2.5) 

^ n =i k ^\ z ) 



where Xj = (xj,yj) and k n (z) := \J z~v? . The complex energy here belongs to 
the resolvent set of H , z G d7 \ [l,oo) ; the function G (-,S 2; ) is defined and 
smooth everywhere outside the hyperplane x\ — , but the sum on the rhs may 
not converge absolutely if the longitudinal coordinates coincide, X\ — x% . Moreover, 
the rhs makes sense also for all non-integer z > 1 , where it gives the boundary value 
of the kernel at the cut; one has to choose properly the branch of the square root 
in k n (z) . We note for further purposes that for real z below the lowest threshold, 
z < 1 , the kernel is strictly positive, Gq(xi,X2',z) > for all mutually different 
X!,x 2 en — cf. [0, App. to Sec.XIII.12]. 

The sought kernel of the full resolvent can be obtained by means of the Krein 
formula App. A], which yields the Ansatz 

(H(a, a) — z)~ l (xi 1 x 2 ) = G (x*i, x 2 ; z) + X(a, a; z)G (xi, a; z)G (a, x 2 ] z) . (2.6) 
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To determine the function A we use the fact that, by definition, (H(a, a) — z)~ l 
maps into the domain of H(a, a) . Writing ip := (H(a, (Tj—z)' 1 ^ , ip := (H — z) -1 y> 
for an arbitrary ip G L 2 (Q) , we get from the last relation 

iX °° e i k n( z )\x~a\ 

ip(x,y) = ip (x,y) H ^ — - sin(ny) sin(n6) ip (a, b) . (2.7) 

71 n =l K n[Z) 

Since ip G D(H ) is smooth at x = a, the generalized boundary values can be 
written as Lj(ip, a) = £j(a)ip (a) . The singularity at the "diagonal" of the resolvent 
kernel for planar regions with a smooth boundary is well known |29| , Chaps. 11 and 
14]. It is not difficult, however, to find it directly: we have 

Cn(a) = lim - — — — — - < 1 + — — sin(ny) sin(n6) > 

\s-a\^o \n\x-a\ \ 7T ^ k n v y; v ;j 

^ ^00 e ik n u ^ ^ 00 ^-(n-z/2n)u 

= — lim — — sin 2 (nb) = — lim sin 2 (nb) , 

7r u^o \nu 1 k n 7r u^o \nu 1 n 

n=L n n=l 

because the differences k~ l — (in) -1 and e tk ^+ n - z / 2n — 1 are Q(n~ 3 ) . The limit is 
easily computed giving £o(a) = —X/2n; this in turn yields 

iX ~ / sin 2 (n6) 1 \ 



Using now the boundary conditions ( |2.2j ) we find that 

A(a,a;z) = A(a,a;z) _1 (2.8) 

with 

w ^ x 1 ^ /sin 2 (n6) 1 \ _ . , . 

A( Q ,a; Z ) = a + - £ ^ - - j =: a - £(«,,) ; (2.9) 

for the sake of brevity we shall sometimes drop a part or all of the arguments in the 
following. 



2.3 Scaling behavior 

Although the modification of these results for a general strip is simple, it is worthy of 
a brief comment. The logarithmic factor in the definition of the regularized boundary 
values refers to a fixed length scale. If we choose the natural scale of our problem 
for it, we replace the factor by In \u(x — a)\ , where v := n/d . In that case £o(d) 
does not change, while C\{d) acquires in addition to the scaling also the additive 
factor In f^J . The formula (|2.8j ) reads then A(a, a; z) = (a — £d(a, z))~ x with 

i ^ { sin 2 (nbv) d 



J + ^GO' (2 - 10) 
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where k n (z) := \J z — (nn/d) 2 . Hence changing the strip width is equivalent to a 

shift in the coupling constant for a fixed d ; one has to replace a by a - ^ In ( . 
In particular, making the strip thin with a preserved corresponds to choosing a 
large positive a in the rescaled problem. 



2.4 The resolvent a la Zorbas 



The resolvent kernel can be obtained even without the boundary conditions Q2.2|) 
if we employ an integral formula due to Zorbas |]32| , Thm.4.1]. According to it, the 
operator family in question can be parametrized by a real 6 ; the function A (a, a; z) 
is replaced by 

A(9, a; z) = (l+e l9 ) _1 — z ) J G (x,a; z)G (x,a;i) dx 



e %6 {i + z) \ G (x,a; z)G (x,a; —i) dx > . 



Using the formula ( |2.5|) together with the identities i — z = k n {i) 2 — k n {z) 2 , i + z 
k n {z) 2 — k n {—%) 2 the rhs can be easily computed; we get 



A(0, a; z) 



-| oo 

E 



Til 



n=l 



_L . Jm ^sing Re / 1 

k n (z) \k n (i) J 1 + COS6 1 \k n (i) 



sin 2 (n6) . 
(2.11) 

Since Re (k n (z) x ) as well as k n (z) — k n (i) are 0{n 3 ) , the series converges and 
can be splitted into the ^-independent and ^-dependent parts. Comparing the two 
parametrizations, we find 



a{6) = F(b) + Sm ° n G(b) 



l + cos^ 



(2.12) 



where 



■y oo 



m ■= - e 



7T 



n=l 



V 



* — — sin (no) — — 

\ 2(n 4 +l) v ; 2n 



G{b) 



-i oo 

E 

7T ^ 



Vn 4 + 1 -n 2 



sin 2 (n6) ; 



the function a clearly maps (— 7r,7r) bijectively onto iR. 



2.5 Spectral properties 

As a rank-one perturbation in the resolvent, the point interaction does not change 
the essential spectrum, i.e., we have a ess (H(a,a)) = [1, oo) with the multiplicity 
2[JE/n 2 ] at a point E , where [■] denotes the integer part. 
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Figure 1: The function £ for different values of b. 



To determine the point spectrum, one has to find poles of the resolvent, i.e., to 
solve the equation 

£(a,z) = a (2.13) 

for z e M. First consider the case z < 1 . It is useful to introduce K n (z) := 
—ik n (z) = y/n 2 — z ; then the function £(a, •) can be expressed as 



£(a, z) 



1 



1 ~ / sin 2 (nb) 

7T ^ \ k„(z) 2n . 



K n (z) 
i oo 

£(a,/) + - £ 



K n (z) 



sin 2 (n6) 



The sum at the right hand side converges because the coefficients at sin 2 (nfe) decay 
like 0(n~ 3 ) as n — > oo . In particular, one can compute easily 



1 

2^ 



ln(2 sin 6) 



(2.14) 



this quantity is positive close to the center of the strip, b 6 , and nonpositive 

otherwise. One has 

d£ 1 ~ sin 2 (n6) 

22 — ^ > °' 



dz 



2tt — i ^(z] 



so the function is monotonously growing in (— oo, 1) ; it diverges at both endpoints. 
We have 

sin 2 b 

£(a,z) = —== + 0(1) as z^l-, 

while on the negative side the definition formula shows that £(a, •) diverges loga- 
rithmically. A natural conjecture about its asymptotics is 



47T 



£(a, z) = - — In - - + — *(1) + 0(1) as 



2tt 



— oo . 



where ^(1) = —7 = 0.57721 (compare with Sec. 1. 5]). Figure 1 suggests that it 
is indeed the case. The asymptotics is independent of b but not uniformly in a as 
the relation ( |2.14| ) shows. 

Another property of £(a, •) is its monotonicity across a halfstrip, 



£(a,z) > ${a!,z) 



iff 



6-^ 
2 



< 



2 



To prove it, we employ the identity sin 2 (n6) — sm 2 (n&') = sin(n(6 + 6')) sm(n(b — b')) 
it shows that 

£(a, z) - £(a', z) = G (0, b + b';0,b- b'; z) , 



for < b' < b < f ; it shows that 
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Figure 2: Eigenvalue plots for different values of a, b. 



Figure 3: Eigenf unctions for different values of b . 



so the result follows from the above mentioned positivity of the free-resolvent kernel. 

Properties of solutions to Eq. ( [2.13 ) follow easily from this discussion: to any 
coupling constant a G M there is just one eigenvalue e(az,a) in (— oo, 1] . The 
function e(-,a) is monotonously growing, 



e(a, a) > e(a , a) iff 
and has the following asymptotic behavior, 

sin o \ 



a > a , 



e(a, a) 



1 - 



ira 



_4 e 2(*(l)-27ra) 



O (et- 



as 



as 



00 



a 



-00 . 



(2.15) 



(2.16) 



Hence H(a, a) preserves the property of the point interaction in the plane, namely 
that it has always a bound state except in the free case; also the asymptotics for 
a — > — 00 is the same. On the other hand, proximity of the boundary pushes the 
bound-state energy up, 



€{a,a) < e{a,a 



iff 



b- n - 
2 



< 



2 



(2.17) 



The resolvent formula ( |2.6j ) provides us also with the (non-normalized) wavefunction 
of the bound state through the residue at the pole, a, a) = G (-, a; e(a, a)) , so 
we have 



ijj(x; a, a) = — 



j co e _ Kn ( e )| x -a| 



7T 



n=l 



we 



sin(ny) sin(n6) . 



(2.18) 



If a decreases, ^ becomes well localized and resembles the Hankel eigenfunction 
of the point interaction in the plane, while in the limit a — > 00 we have 



sino 7r 



- |aj— a| \/n 2 — 1 



n=2 



rr 



sin(ny) sin(nfe) 



the leading term here is the product of Xi(y) with the eigenfuction of the one- 
dimensional attractive point interaction of the strength — (2/na) sin 2 b ; recall that 
in view of the preceding section this corresponds by rescaling to the limit d — *■ . 
These results are illustrated on Figs. 2 and 3. 

There are no eigenvalues embedded in the continuous spectrum. The function 
£(a, •) is defined everywhere in (l,oo) except at the thresholds, z = m 2 , but it 
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is now complex because K n (z) is negative imaginary for the open channels. In 
particular, 

, . 1 sin 2 (n6) 

lm£(3,z) = - E -TTT- > ( 2 ' 19 ) 

71 n =l K n\Z) 

for a non-integer z > 1 , so the function a— £(a, •) has no zeros there. The argument 
does not apply if z — m 2 because the corresponding k m (z) = and none of the 
expressions at the rhs of (|2.6|) makes sense. However, one can check directly that 
the singularities cancel mutually. Around z = m 2 , the free resolvent kernel behaves 
as 

n,-. - \ i sm(myi)sm(my 2 ) 

G(xx,x 2 ;z) = Err(z), 

7T \J z — m z 



where Err(^) stands as a shorthand for 1 + 0(yz — m 2 ) , and 

i sin 2 (m6) 



/ — * \ ~ 0111 \ 7 //-\ / I 

A(a, a; z) = a — — / = + U(y z — m 



71 V 



z — m 



with a := a + tt— — 1 Y,. n ^ m ( rV"^ — k^- V Hence the full-resolvent kernel 
behaves asymptotically as 

ia sin(m?/i) sin(m?/2) ( ~ i sm 2 (mb) ^ 1 ^ ^ 



7r \J z — m 2 \ vr \fz 



m 



and cannot have a pole at z = m 2 . In conclusion, no real z > 1 is an eigenvalue 
of H(a,a) . On the other hand, the resolvent has resonance zeros in the complex 
plane as we shall show in Section 2.7 below. 



2.6 Scattering 

The existence of the wave operators is easy to establish because the point interaction 
is a rank-one perturbation in the resolvent, and therefore Birman-Kuroda theorem 
applies [^, Sec. XI. 3]. What is more interesting, however, is how the scattering 
matrix looks like. To find it we employ the relation ( J2.7|) in combination with ( |2.8| ) 
changing slightly the notation: to any ip G D(H(a,a)) and a nonreal z there is a 
unique decomposition 

i/>(x) = ^z{x) H G (x,a;z)ip z (a) 

a - £{a,z) 

with ip z e D(H ) and (H(a, a) — z)ip = (Ho — z)ip z . If we choose 

i/f g := e lkniz)x - £x2 Xn {y) (2.20) 
for ?/> 2 , the corresponding ip =: i\f belongs to D(H(a,a)) for all e > and 
{{H(a,a) - 2)^ £ )(x) = 2e(2 £ x 2 - 1 - 2ik n {z)x)ip E z (x) . 



10 



Since the rhs converges in L 2 sense as z approaches the real line, ip e G D(H(a, a)) 
for z G [1, oo) and the last relation holds again. Of course, ip £ ceases to be square 
integrable as e — > 0+ , but the pointwise limit exists and 

p ik„(z)a 

^(f) = e ik ^ x xn(y) + ~cT~* V G ^ z ^(b) ■ ( 2 - 21 ) 

a — £(a, z) 

This function is locally square integrable, satisfies the appropriate boundary condi- 
tions and solves (H(a, a) — z)ip = as a differential equation. Substituting for G 
we obtain easily the reflection and transmission amplitudes, r nm (z) and t nm (z) , 
respectively, for the case when the incident wave corresponds to the n— th trans- 
verse mode and the scattered one appears in the m-th mode; they are given by the 
relation 

(+ (?)-X \J(km-k n )a _ , x -i(k n +k m )a _ 1 sin(w6) sin(m&) 

7T fc m (z)(a - 4(a, z)) 

It is necessary to assume, of course, that both the involved channels are open, i.e., 
z > max(n 2 , m?) . 

Since the point interaction gives rise to a rank-one perturbation in the resol- 
vent, the Kato-Birman theory yields also the completeness of the wave operators, 
i.e., unitarity of the S-matrix. In particular, the "diagonal" part of the unitarity 
condition can be expressed in terms of the transmission and reflection coefficients 

as 

^ ] /c m (|t nm | + |r nm | ) k n , (2.23) 

m=l 

and has an obvious meaning of conservation of probability flow (taking into account 
different particle velocities in different channels). It is not difficult in the present 
case to check the identity ( |2.23| ) directly. Since t nn = 1 + r nn e~ 2tkna , it reduces to 

k n Re \ T nn e J + 2 ^ ] k m \r nrn \ , 

m=l 

substituting from ( |2.22| ) and multiplying both sides by \a — £(a, z)\ 2 , we arrive at 
the relation 

sin 2( m6 ) 

Imf(0 ' 2) ^£^r- 

which is nothing else than (|2.19|) ; recall that the closed channels do not contribute 
to the imaginary part of £ (a, z) . In the same way one can check the remaining parts 
of the unitarity condition which read 

^ ] km(tnmtsm TnmTsm) $nskn 
m=l 
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for n ^ s and 

Eh It r p2ik m a , T -2ik m a\ _ n 

m=l 

since the "elementary block" of the S-matrix describing the transitions between the 
n-th and m— th channel is of the form 



,r __ km ( trim ^nm \ .) i 

Onm — \l -j— I ~ 7 ) , ■ ' 

n^n \ ' nm b r, 



where the tilded quantities in the second row are obtained by changing the sign of 

an p rp • rp p 2 2 {kjl H~&m pj- p 

2.7 Resonances 



Every square root in the definition of the channel momentum k n (z) := V z — n 2 
gives rise to a cut; hence the free resolvent kernel ( |2.5|) as well as other quantities 
derived from it are in general multivalued with infinitely sheeted Riemann surfaces. 
In particular, the pole condition (|2.13|) may have solutions on the other (nonphysical) 
sheets which produce resonances. 



For simplicity we introduce q n (z) '■— V z — n 2 which takes values always in the 
upper complex halfplane or at the positive real halfline. On the (N + l)-th sheet 
we have then 

f -q n (z) ... n < N 

\ q n {z) ... n > N 



kjz) 



Denoting by ^jv(o> •) the corresponding branch of the function £(a, •) we can write 
the pole condition (|2.13|) explicitly as 



a — Re^Ar(a, z) = a 



+ - V ( -Im (—^—\ sin 2 (nb) + — 



2n. 



+ - E ( Im ( ) sin 2 (n6) + ^ ) = 



(2.25) 



1 JI / 1 



- Im^a, z) = - V Re —— sm 2 (nb) 

-- V Re sin 2 (n6) = 0. 

For a weak coupling, there is generically one resonance pole close to each threshold, 
with exception of the lowest one; the resonance is absent if (iV+l)6/7r is integer so 
that the incident wave has a node at the impurity and does not feel the perturbation. 
To show this more explicitly, let us rewrite £at(<J, z) = £ N (a, z) + ^ N (a, z) , where the 
two terms sum contributions from the evanescent and propagating modes, n > N+l 
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and 1 < n < N , respectively. Suppose we switch the latter on with the help of an 
additional parameter, i.e., we look for solutions of the condition 



F(z, V ) := a-^(^z)- V ^(a } z) 



0. 



(2.26) 



If r] = , we can repeat the argument of Sec. 2.5: we find that the function £ N (a, •) 
is monotonously growing between — oo and oo when z runs over the interval 
(— oo, (iV+1) 2 ) ; hence to a given a there is just one z° N (a) such that F(z%(a), 0) = 
. Moreover, the leading behavior of £ (a, •) as z — > (iV+l) 2 — is again given by a 
single term, so the "eigenvalues" behave as 



4(a) = (N+lf 



/sin 2 ((iV+l)fc)' 



0(a 



(2.27) 



in the weak-coupling case, a — > oo . 

For a nonzero r] the condition ( 2.26Q can be solved perturbatively by means of 
the implicit-function theorem. We have 



OF 
dr] 

dF 

dz 



N 



(4H>°) 



-r (4(a)) = - E 

7T 



sin 2 (n6) 



n=l 



'(iv+i; 



n 



= + f + c(0> 

2 2zra ' 



- E 

2tt z ^ 



7r 3 a 3 



■;i=Af+l 



sin 2 (nb) 

?iv(4(«)) 3 = |sin((iV+l)6))| 



+ C?(a°) 



dividing these quantities we obtain —(dz^/d^^Q . Moreover, the remainder term 
coming from d 2 z n N /drj 2 is 0{a~ i ) , so for sufficiently large positive a we may use 
the expansion up to 77 = 1 , obtaining an asymptotic formula for the resonance-pole 
position, 



z N la 



2i 



sin((iV+l)6)| £ 



sin 2 (rT,6) 



n 3 a 3 



n=l 



'(N+l) 



+ — + 0(a 

2in 



(2.28) 

the real shift appearing in (|2.27|) has been absorbed here into the error term. In 
particular, the resonance width behaves in the weak-coupling case as 



T N (a) 



-2ImzAr(a;) = 4 



sin((JV+l)6)| ^ 



sin 2 (nb) 



7r 3 a 3 



^! ^(AT+l)2-n 2 



+ 0(a- A ) . (2.29) 



If the coupling is not weak, resonance poles can be found from a numerical solution 
of the conditions (|2.25|) ; an example is shown on Figure 4. The pole trajectories 
can be followed for all values of a , however, only weak perturbations produce a 
substantial resonance scattering effect because the pole residues decrease rapidly 
with the coupling strength; this also is shown on Figure 4. 
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Figure 4: Pole trajectory in the z-plane. 



3 A finite number of point interactions 

Our next aim is to extend the above results to the situation when there is a finite 
number iV of point interactions in the strip Q . We suppose that their positions 
are dj := (aj, bj) and the coupling constants ctj ; for the sake of brevity we denote 
a := {di, . . . , <2at} and a := {«i, • • • , ajy} . The corresponding Hamiltonian H(a, a) 
is again given by ( |2.3| ) with the boundary conditions ( |2.2p replaced by 

Li(i/j,aj) + 2irajL (i/;,dj) = 0, j = 1, . . . , N . (3.1) 

As in the one-center case, any of these point interactions may be switched off by 
means of the limit atj — > oo . 

3.1 The resolvent 

We shall again start from derivation of the resolvent. Since the deficiency indices of 
the initial symmetric operator are now (iV, N) , a natural Ansatz is 

TV 

(H(a,a)-z)~ 1 (x 1 ,x 2 ) = G a {x u x 2 ]z) + Y ^jkG (x 1 ,a j ;z)Go(d k ,x 2 ;z) , (3.2) 

j,k=i 

where Gq(-, •; z) is the free-resolvent kernel (|2.5|) . Functions ip G D(H(a,a)) and 
■00 ^ D(Ho) introduced as above are then related by 

TV 

ip(x) = ipo{x) + Y ^jkG (x,a j ;z)^ (d k ) (3.3) 
j,k=i 

and one can compute easily the boundary values at the interaction support, 



N \ 
2tt 



L (^,a m ) = - Y i^; 5 jm^o(d k ) 
j,k=i 



Lx(ip,d m ) = -ip (a m ) + Y ^jk -Sj m Y 



' i ~ / sin 2 (nbj) 1 



i 7T 

j,fc=l \ n 



^ \ fc n (;z) 2in 



+ (1 - 5 jm )G (dj, a m ; z) \ipo{d k ) . 



Substituting this into ( |3.1|) and using the fact that ipo(dk) are independent quanti- 
ties, we find 

X(a,a;z) = A(a,a; zy 1 , (3.4) 
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where A (a, a; z) is now an iVxAf matrix given by 

i ~ / sin 2 (nbj) 1 \ 
^ k n {z) 2m J ' 

(3.5) 

2 oo e ife„(z)|oj-a m | 

-G Q (aj, a m ;z) = TTT~\ sin(n6 j ) sin(nfe m ) , j ^ m ; 

^ n=l fc nl^J 

we shall again drop the arguments occassionally if they are clear from the context. 
Notice that the scaling argument of the Section 2.3 can be easily adapted to the 
present situation; it shows that changing d is equivalent to the simultaneous shift 
of all the coupling constants <x, on — i In J in the rescaled strip. 



A 



Aj m 



3.2 The discrete spectrum 

A finite number of point interactions represents a finite-rank perturbation in the 
resolvent, and therefore the conclusion of Sec. 2. 5 about the essential spectrum re- 
mains valid, a ess (H(a, a)) = [1, oo) . The discrete spectrum is again determined by 
poles of the resolvent coming from the coefficients Xjk in ( |3.2| ); they are given by 
the condition 

det A(a, a, z) = . (3.6) 

Comparing with the case AT = 1 , it is now slightly more complicated to determine 
the eigenfunctions. One can adapt the argument from |], Sec.II.l]. Suppose that 
H := H(a, a) satisfies Hp = zip for some z G M. We pick an arbitrary z' G p{H) ; 
then in analogy with ( |3.3| ) there is a vector tp G D(H Q ) which allows us to express 
the eigenvector p as 



N 



(3.7) 



ipo + ^2 d j G o(-,aj;z') , 

where the coefficients are given by dj := Ylk=i(^( z ') ~ 1 )j&'0o(Qfc) and the relations 

(H - z')ip = (H - z')<p = (z - z')<p 
are valid. Applying (H — z'y x to the last identity we obtain 



[z — z 



N 



(H - z 



ipo+ d i( H o ~ z ') lG o(-, z) 

.7 = 1 



and this in turn yields 

(H - z)ip 



N 

y z — z') ^2 djGo(-, 3j] z') 

3=1 



(3.8) 
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If z < 1 , the resolvent (Ho — z) 1 exists and may be applied to both sides of the 
last relation giving 

N 

1>o = XX' (^o (•,%■; z) -Go(-,aj\z')) ; (3.9) 



we have employed here the first resolvent identity. Substituting into ( ft.7| ) we get an 
expression for the (non-normalized) eigenfunction 

N 

<f(x) = djG (x, a,-; z) . (3.10) 
i=i 

To determine the coefficients, we notice that the relation ( |3.9j ) together with ( |3.5| ) 
gives 

N 

V'o(oj) = X ( A ( 2 ')jm - A(z) jm )d m ; 

m=l 

on the other hand, the above mentioned expression of dj is equivalent to ipo{dj) = 
J2m=iM z ')jmd m - Hence 

N 

A{z) jm d m = , (3.11) 

m=l 

i.e., d : = (d\, . . . , aV) is an eigenvector of A(a, a, z) corresponding to zero eigen- 
value. Inverting the argument as in one can check that any solution to ( [3.11 ) 
determines an eigenvector of H(a, a) . 

The next question concerns the existence of solutions to Eqs (|3.6j ) and ( 3.11 



For convenience we introduce again K n = —ik n (z) = \fn 2 — z which is positive for 
z < 1 . We have then 

1 ~ f sin 2 (nbj) 1 \ 



Aii = % - - £ 



^ „=i 



« n (z) 2n 



and an analogous expression for the nondiagonal elements. If z — > — oo the matrix 
behaves as 

A(a,a,z) = i- ln(-|) I + , 

so all its eigenvalues tend to +oo as z — > — oo . On the other hand, for z — > 1 — 
we have 

A(a,a,z) = 7=— ^i + 0(1), 

where Mi := (sinfej sinfe^j^^. This matrix has, in particular, an eigenvector 
(sin fei,..., sin 6jv) corresponding to the positive eigenvalue J2j=i sin 2 > an d there- 
fore one of the eigenvalues of A (a, a, z) tends to — oo as z — > 1— . The elements of 
A (a, a, z) are continuous functions of z , hence the same is true for its eigenvalues. 
It follows that there is an eigenvalue which crosses zero, i.e., H(a,a) has at least 
one eigenvalue. 
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Figure 5: Eigenfunctions for three point perturbations, (a) The ground state, (b) 
The second excited state. 



In fact a stronger claim can be made concerning the eigenvalues of A (a, a, z) 
By a straightforward differentiation we find 



d l 00 e -\aj-a m \Vn 2 -z 

— A(z) jm = - — ^ n 2 _ ^3/2 V 1 + l a J _ a mWn 2 - zj sm(nbj) sm(nb m ) . 

The matrix function A(-) is monotonous if for any c G <U N the quantity 
— (c, A(z)c) = ~ — J2(n 2 - z)~ 3/2 Cj sm(nbj) c m sin(nb m ) 

dZ Z71 n=l j,m=l 

x e -K-a m |^?^ ^ + ^ _ am |Vn 2 -z) 

has a definite sign (is nonpositive in our case). The expression on the right hand 
side tells us that this is true provided the function / : f(x) = e~ K ^(l + k\x\) is of 
positive type for any k > . The last named property follows from the identity 

2 k 3 f e ipx 
(l + K |x|)e~ KN = / dp, 



because by Bochner's theorem Sec. IX. 2] a function is of positive type iff its 
Fourier image is positive. In fact, since the measure in the last integral is point- 
wise positive, 4- A(z) is even strictly positive; it means that all the eigenvalues of 
A (a, a; z) are decreasing functions of z . 

The most important consequence of this fact is the nondegeneracy of the ground 
state. To prove it one has to check that the lowest eigenvalue of A(z) is simple 
for any z G (—00, 1) , which is equivalent to the claim that the matrix semigroup 
I e -*A(z) . t > 0} is positivity preserving |25|, Sec. XIII. 12]. The last property is 
ensured if all the nondiagonal elements of A (z) are negative; we have A(z)j m = 
—G (dj, a m ; z) by ( |3.5D so the desired result follows from the positivity of the free- 
resolvent kernel. The coefficients may be therefore chosen of the same sign for the 
ground state; in fact, as positive because dj = would mean that the eigenfunction 
( |3.10|) is smooth at t J QjJq SO ;he corresponding interaction is absent, <x, = 00 . 



3.3 Embedded eigenvalues 

Let us now ask whether the continuous spectrum of H := H(a,a) may contain 
eigenvalues if > 1 . Suppose that Hip = zip for some z > 1 . We employ 
again the expression ( |3.7| ) for the eigenvector and write ip as a series, ipo(x) = 
12™=i 9n{ x )Xn(y) ) with the coefficient functions g n G L 2 (M) . Substituting this into 
( |3.£| ) and using the fact that {Xn} is an orthonormal basis in L 2 (0,d) , we obtain 
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the following system of equations, 



„• N ik n {z')\x-a,j\ 



-gjx) - K(z) g n (x) = ^(z-z)}2 djXn(bj; fe ^ 

for n = 1,2,... ; the Fourier-Plancherel operator transforms it to 

(p 2 -z + n 2 )g n ( P ) = —j: d ^) p2 _ zi + n2 . (3.12) 

If g n G L 2 the same is true for g n ; this is impossible if 

z > n 2 (3.13) 



and the rhs of ( ft.12;) is nonzero at ±p n , where p n : = \Jn 2 — z , since g 2 would have 
then a nonintegrable singularity. It is clear that it is the factor p 2 — z + n 2 which 
matters, because z' belongs to p(H) by assumption. 

Hence the rhs of ( P-12 ) has to be zero; we want to conclude that g n = . If 



N > 1 and ctj are not the same, it might happen that the rhs is not zero identically. 
However, we can choose a common phase factor to be put in front of the sum; then 
the square integrability requires 



N 

E 



for an arbitrary a . If all the a,j are mutually different (mod 2irp~ l ) it follows 
that djXn{bj) = for each j . On the other hand, if some of them coincide we find 
J2j djXn(bj) = where the index runs through the values with the same longitudinal 
coordinate aj , and therefore g n = again, i.e., g n may be nonzero at most if some 
Oj differ by multiples of 27rp~ 1 . 

Consider now an arbitrary g G L 2 (M) and n satisfying the condition ( 3.13|) . 
Using (|3.7|) and (|2.5|) we find 



i 



2k n (z j .. | 



where the inner product on the rhs refers to L 2 (M) and g n in the first term can 
be expressed by Q3.12 ). If djXn(bj) = for each j , the right hand side is zero. In 



the exceptional case mentioned above we use the fact that the Ihs is independent of 
z' . The explicit expression for dj together with the asymptotic behavior of A show 
that dj — > as z' — ► — oo ; the same is true for the inner product in the second 
term as well as for (g,g n ) ; together we find (gXn,<f) = again. We conclude that 
z > 1 cannot be an eigenvalue corresponding to an eigenvector from the subspace 

® [ £[L 2 {IR)®{x n }. 
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Figure 6: An embedded eigenvalue due to symmetry. 



On the other hand, the condition ( |3.13| ) in the above argument is crucial; in 
the case N > 1 the operator H(a, a) can have embedded eigenvalues with eigen- 
functions in the orthogonal complement of the mentioned subspace. As the simplest 
example, consider a pair of point perturbations with the same coupling constant 
a placed at a\ := (0,6) and a\ := (0, 7T — 6). The eigenvalue problem can be 
divided into the part symmetric and antisymmetric with respect to the strip axis, 
the antisymmetric part being obtained by scaling of the single-center problem with 
a := (0, 2b) and coupling constant a— In 2 . The scaled eigenvalue tends to 4 as 
a — > oo , hence it is embedded in the continuum for all a large enough. This is illus- 
trated on Figure 6. In the same way, one can construct other examples of embedded 
eigenvalues. Their common feature is the existence of a symmetry which prevents 
the (energetically allowed) decay of the eigenstate; a violation of the symmetry 
turns these eigenvalues into resonances. Recall also that embedded eigenvalues due 
to symmetric obstacles in Neumann waveguides have been treated recently in [|TT|. 



3.4 The limits of strong and weak coupling 

If all the point interactions under consideration are strong, i.e., having the a/s 
large negative, one naturally expects the corresponding bound states to be strongly 
localized and weakly influenced both by the other perturbations and by the bound- 
ary. To show that this is indeed the case, let us write the matrix A(a, a; z) in the 
form 



A(z) 



1 , ( z 
atj + — In - - 

3 4ir I 4 



I + 



1 , ( z 
atj + — In — - 

47T V 4 



5j k 



-i 



A(z) 



where A(z) is the remainder matrix, which is independent of a and has a bounded 
norm as z — > — oo . Given a finite energy interval I C (— oo, 1) one can always 
choose them's large enough negative so that the two matrices in the above product 
are regular, and no eigenvalues are contained in / . Consequently, the roots of 
Eq.( |2.13D in the strong-coupling limit are situated in the region where A(z) is 
dominated by the diagonal part; then there are exactly iV eigenvalues, including a 
possible degeneracy, and 



€j{a,a) ~ 



4e" 



-47TO, 



as 



max a,- 
i<j<N J 



-oo . 



(3.14) 



One can conjecture that the asymptotics is in fact the same as in fl2.16|) . The 
corresponding eigenf unctions coincide in the leading order with the strongly coupled 
one-center eigenf unctions of Section 2.5. 

For a weak coupling the situation is different. We again restrict our attention 
to the case when all the interactions are weak, i.e., the corresponding a/s are large 
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positive. We denote A := diag(«i, . . . , ary) and use the decomposition 



A(z) 



A - 



7Ty/ 1 — Z 



Mi 



A - 




TTy/1 — Z 



where A(z) is a remainder independent of a , whose norm is bounded as z — > 1— . 
In the same way as above, to a given z < 1 one can always choose the o^'s suffi- 
ciently large so that the two matrices are regular and no eigenvalues are contained 
in (—00, zq] . In distinction to the previous case, however, A(z) is asymptotically 
not a multiple of the unit matrix but rather a rank-one operator on as we have 
seen in the previous section. Hence only one eigenvalue of H(a, a) can approach 
the threshold of the continuous spectrum; in combination with the above argument 
this means that if all the point interactions are weak enough, H(a, a) has a single 
bound state. 

This is analogous to the behavior of interactions whose support has a nonzero 
measure. Recall that a two-dimensional Schrodinger operator with a weak potential 
V has always one bound state provided V is attractive in the mean and satisfies 
some decay restrictions fl26| |, and the same is true in a tube, whatever dimension its 
cross section is, provided the potential is replaced by its projection on the lowest 
transverse mode ||I2"| . 

To find the corresponding bound state in our case, we have to solve the spectral 
problem for Mi + t]A , where i] := — ity/l 
that approaches J2f=\ srn2 bj as z — ► 1—. 
yields the expression 



- z ; we are interested in the eigenvalue 
An elementary perturbative argument 



N 

3=1 



sin 2 bi 



7T 



Y/. 



N 



sin 



2 h 



+ 0( V 2 



putting it equal to zero, we obtain the energy of the weakly bound state, 



e(a, a) 



fetish^ 

EN 2 l 

• =1 OLj sin bj 



(3.15) 



as mm.i<j<N a j 00 • Since the range of Mi is spanned by (sin b%, . . . , sin &ry) , 
the corresponding asymptotic expression for the eigenfunction is 



ip{x\ a, a) 



smy 



■ 2 u N 



(Etisin 2 6 fc ) 



^ e -v^|x- % l sin 2^ 



00 N g— y/n 2 — l\x— Oj| 

+ ^ sin(n?/) ^ y = — sin(n6j) sinfoj . 



n=2 j=i v 7 ^ 2 - 1 



The leading term is again a product of Xi{y) with a linear combination of the 
eigenfunctions of one-dimensional point interactions placed at a,-, j — 1, . . . , N . 
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3.5 More about bound states 



The spectral condition ( |3.6|) makes it possible to derive also other properties of the 
discrete spectrum for the A-center Hamiltonian. We have seen that it has iV 
bound states and just one bound state in the strong- and weak-coupling limits, 
respectively. One may ask, more generally, what are the regions in the space M N of 
coupling constants, where the dimension of the discrete spectrum is n = 1,2, . . . , N . 

First consider the two-center case. As we already mentioned detA(a,a;z) tends 
to +00 as z — * —00 . On the other hand, the behavior of this quantity as z — > 1 — 
depends on the parameters a.j, Sj, j — 1, 2 ; the limit may be either +00 or —00 , 
i.e., the equation det A(a, a; z) = can have one or two roots. The corresponding 
areas in the (ax, a 2 )-plane are divided by the line ^2(^1, ^2) = , where 

p ( \ • 2l -2, sin 2 61+ sin 2 6 2 2 . 2 , . 2 , 

r2\a\,a.%) := ol\ sin o 2 + a 2 sin b\ H — \a,\— a 2 | sin o x sm o 2 ; 

27T 7T 

there is only one bound state if Pi(a.2, olq) > and two in the opposite case. 
Inspecting the expression more closely we see that for a fixed a\ it represents a 
family of lines with the common intersection 

4|ai— 02] sin 2 b\ — 1 1 
_ 2^ ' _ 2tt_ 

Obviously, we can change the number of bound states for fixed coupling constants 
by changing the slope of the line (it depends on bi,b 2 only) or shifting the line in 
the «i-direction changing \a\ — a 2 \ . 

Similar conclusions can be made for a larger number of perturbations. Although 
parametric expressions of the hypersurfaces in M N dividing regions with different 
number of bound states are easy to derive investigating \im z ^i_ det A(a, a; z) , the 
formulas are somewhat lengthy and we do not write them down here. We restrict 
ourselves to the observation that 

P/v(«i, . . . , a N ) = Pjv-i(«i, • • • , a N -i) a N + Qn-i(chi, a N -i) 

for some polynomial Qn-i, which means that the dividing hypersurface for N 
perturbations is singular as a function of ai, . . . , ajv-i at those points for which 
Pat_i («!,..., ctjv-i) — 0, i.e., on the previous surface. Number of bound states 
is changed by one each time when the surface in the coupling-constant space is 
crossed. It means that having found the number of bound states for a given setting 
of N—l centers, we can decide whether adding another center increases the number 
of bound states or not: it does if P^(a\, . . . , a^) < , while in the opposite case it 
does not. 

Another question concerns possible degeneracy of the discrete spectrum. For 
the sake of brevity, we write the matrix ( |3.5| ) as 

A(a,a;z) jm = 6 jm (atj- £j{z)) + (l-5 jm )g jm (z) 
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with £j(z) := £,(dj,z) and gj m (z) := —Go(aj,a m ;z). Since —gj m ( z ) > 0, the 
maximum degeneracy is N—l ; in particular, the discrete spectrum is always simple 
for N = 2 . This is not surprising; we know from Section 3.2 that the ground state 
is nondegenerate. On the other hand, degeneracies may occur. Consider the case 
N = 3 with 01,3 := f±a, ~) and a 2 '■= (0, b) , and fix an energy z < 1 . We 
have (712(2:) = tot^) for any b G (0, tt) . If 6 = | this value is by Q2.5| ) strictly 
greater than #13(2) ; on the other hand, lim^ gniz) = , so there is a b G (O, |) 
for which all the three gj m {z) have the same value. Choosing now the coupling 
constants a>j in such a way that ctj — £(<%; z) = gj m (z) , j = 1, 2, 3 , we find that z 
is an eigenvalue of multiplicity two. 

It should be noted, however, that degenerated eigenvalues occur rather excep- 
tionally. The matrix A(a, a; z) has 3iV real parameters, because one of the coordi- 
nates dj may be chosen arbitrarily, while the number of different 2x2 submatrices 
is (^J 1 ) • 111 particular, the number of the available parameters is always less than 
the number of matrices to be annulated for N > 6 . 

Suppose now that we fix N—l coupling constants and let the remainig one run 
from —00 to 00 . The whole discrete spectrum moves up at that but in a peculiar 
way; if we do not hit a degeneracy point by a chance, we observe avoided cross- 
ings between the energy levels of the (iV — l)-center Hamiltonian and the running 
eigenvalue corresponding to the chosen perturbation. If all the fixed <x,- represent 
a strong coupling, so that they meet the graph of the chosen level at large negative 
energies where it is steep, and if the crossings are narrowly avoided, we reproduce 
in the present situation the cascading phenomenon studied in fl7 |. 

To explain the mechanism responsible for this behavior, consider the case of two 
perturbations with coupling constants a := a>\ and a 2 . If one of the perturbations 
is absent, the other one has a single eigenvalue, ei(a) and e 2 (a) = e 2 , respec- 
tively. Now we let a running and look how the eigenvalues influence each other, in 
particular, when \a\ is large. One of the solutions to the condition (|3.6|) is 



where 



e 2 (a) = e 2 + + 0(a~ 2 ) , (3.16) 

cto— a 

<H{z) .= -g^- > 0, = ^ " £(e 2 ) < 

the first term on the rhs of the last expression corresponds to the crossing point, 
ei(a) = e 2 , and the second one is the shift due to level interaction. In a similar 
way, the running eigenvalue changes to 

e / x _ e , x + ffi2(ei(a)) 2 

1 ^i(ei(a))[^ 2 (ei(a)) - a 2 \ - 2(g 12 g' 12 )(e 1 (a)) ' 

Using asymptotic properties of ^ together with (|2.16|) , we find that 

4-7T 

ei (a) w e 1 (a) + — (gn^a))) 2 e W)-*«) (3.17) 
a 
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Figure 7: The cascading effect for a cluster of eigenvalues 



as a — > — oo ; recall that g\i{e.\ (■)) itself is exponentially decaying. On the other 
hand, there is no reasonable asymptotics for a — > oo because the cross terms #12(2) 
can no longer be regarded as a perturbation there; in fact, the perturbed eigenvalue 
may disappear in the continuum. 

Nevertheless, if 2gi 2 g' 12 remains small in a wide range of energies, there is an 
interval (a', a") above the crossing point where the first term in the denominator 
dominates and we have 

f \ , N cAeAa)) 
eAa) « eAa) + 



6(ei(a)) - "o 

with a := a2 + (2g 1 2g' 12 /^i)( e i( a )) j the correction term is positive. We see therefore 
that away of the crossing region the two eigenvalues follow closely the two branches 
of the "decoupled" spectrum graph. It is also clear that the smaller are the terms 
gi2 in the above formulas the closer come the two curves together; since the function 
Go(ai,a2; z) at a fixed energy decrease exponentially with the distance of the two 
points, the most profound cascading effect may be expected when the impurities are 
far apart. Moreover, if the perturbations produce a multiple eigenvalue or a cluster of 
almost identical simple eigenvalues, the cascading means that one eigenvalue leaves 
the cluster and one joins it; this is illustrated on Figure 7. 

3.6 Scattering 

The existence and completeness of wave operators follow again from the Kato- 
Birman theory. The argument of Section 2.6 can be easily modified: the generalized 
eigenf unction ( |2.21|) at a non-integer z > 1 has to be replaced by 



= e lk ^ x Xn(y) + E (A(z)- 1 ) lk G (x,a 1 ;z)e^ a *Xn(h) (3.18) 
j,k=i 

with the incident wave in the n-th channel. Consequently, the reflection and trans- 
mission amplitudes from the n-th to the m-th channel are 



? ' n.m. \Z) 



tnm (^) 



i (^A( s ^ n ( m ^j) sin(n6 



71 jj^=L J k m (z) 

(3.19) 

7T j k=1 Km[ z ) 



The unitarity condition now reads 



m=l 
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Figure 8: The conductance plot for a pair of point perturbations. 



(3.20) 

m=l 



because S nm is given again by ( |2.24j ), where the tilded quantities are obtained by 
mirror transformation, i.e., by changing each perturbation longitudinal coordinate 

dj tO — CLj . 

The knowledge of the reflection and transmission coefficients makes it possible 
to express quantities of a direct physical interest, in the first place the conductance 
which is given by the Landauer formula 

g{z) = ^ E t^IW*)! 2 , ( 3 - 21 ) 



h n,m=l 



where t nm (z) are the coefficients ( 3.19|) . An example is shown on Figure 8; we see 



that the perturbations deform the ideal steplike shape. The most remarkable feature 
are the sharp resonance peaks and dips which approach the channel thresholds in 
the weak-coupling limits. The corresponding resonance poles can be found as in 
Section 2.7 but we shall not go into the details. 

The fact that our model is explicitly solvable makes it a useful tool to study 
various methods in which the scattering matrix for a system of perturbations is 
constructed through a factorization from single-obstacle objects. In practice such 
techniques involve a restriction on the number of evanescent modes taken into ac- 
count between subsequent scatterers. A detailed discussion of this situation in our 



framework is presented in another paper \W\; it shows that a particular caution is 
required in the resonance regions. 



4 A cylindrical strip in axial magnetic field 

Before proceeding further, let us show how the above reasoning can be modified for 
a different geometry: we shall suppose that the strip is coiled into the form of a 
cylinder and placed into a homogeneous magnetic field parallel to the cylinder axis. 
Here we consider the situation when the longitudinal coordinate is locked cyclically 
so the resulting cylinder is of a finite height; the case of an infinite cylinder which 
was discussed recently on a heuristic level in will be treated in the next section. 



For a finite-height cylinder the configuration space is compact and may be 
mapped onto the rectangle 

tt := [0,2tt) x [0,d] 
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with the cylinder perimeter equal to 2ir on a suitable length scale. The vector 
potential at the cylinder surface is A = (|£>,0J , so 



cf) ■= n B = 2txA (4.1) 
is the magnetic flux through the cylinder. 

4.1 The free Hamiltonian 

The free problem can again be solved by separation of variables. It is convenient 
to include the magnetic field into boundary conditions: the unperturbed operator 
H ((p) acts then as Laplacian on Q satisfying the condition ( [2.1| ) for all x G [0, 2ir) 
together with 

^(0+,y) = e-^arr- y), ^(0+,Z/) = e^^(2n-,y) (4.2) 

for all y G [0, d] . Its spectrum is purely discrete (i.e., finitely degenerated) and 
consists of the eigenvalues 

e mn := (m + A) 2 + (un) 2 , m G 2Z , n = 1, 2, . . . ; v := ^ . (4.3) 

Lb 

In particular, if v 2 is irrational the spectrum is simple for 2A G" /Z and twice 
degenerated otherwise (with the exception of the ground state for A G 2Z ) . The 
corresponding eigenfunctions are of the form rj m ® Xn > where \n are the transverse- 
mode functions used above and 

rj m (x) := -Le^+^s, mE Z. 

It is straightforward to write the unperturbed resolvent kernel 

^ _ _ , , 1 ~ =2, e i ( m+A ^ :Cl - a;2 )sin(z/ra?/ 1 )sin(z/ra?/2) , v 

G o (aci,a? 2 ;0,2;) = -j L 2^ 1 ■ ( 4 - 4 ) 

™ n=l m=-oo e mn 2 

The rhs has to be handled with caution because the series is not absolutely conver- 
gent. Fortunately, one can get rid of the double summation by evaluating the inner 
series: we have 

oo e imx ^ / e %{A-i^f){-n-x) e i(A+i^y)(TT-x) \ 

m ±L (m + A) 2 + 7 ~ 27^ ^sin(7r(A-i^/7)) ~ ^RI+V^j) J 

for x G [0,27r) p4|, 5.4.3.4], while for x G (— 27T, 0) one has to replace n — x by 
— 7r — x ; this yields 

sinh((27r— \xi — X2I) \/u 2 n 2 — z) + e 2i,?7rj4 sinh(|xi — a^lV^ 2 ^ 2- 2) 



cosh(27rv^ 2 n 2 — 2) — cos(27rA) 

(4.5) 

sin(i/ra/i) sm(i/ny2) 



x 



v 2 n 2 — z 



where r\ := sgn (xi — £2) 
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4.2 The perturbed resolvent 

As before we add now to the operator H ((p) a finite number of point perturbations 
with the coupling constants a\, . . . , placed at a%, . . . , 6 ft ; the corresponding 
Hamiltonian H(a,a,<p) is again determined by the boundary conditions ( |3.1|) . To 
find its resolvent we repeat the above argument writing a general vector from the 
domain of H(a, a, 0) in the form Q3.3|) . In view of ( |4.5|) , we have 

G (x,dj) 1 1 ~ wx(ynyi) sm(pny 2 ) , . 

hm - — : — = — hm - — : — y coth(7ri/n) . 

\x-a,j\-*o m\x — dj\ 7r y-'bj la\y — bj\ ^ n 

Since coth(irun) tends to one exponentially fast as n — > oo , it may be neglected; 
Using further the identity 



sm(unyi) sin(vny2) 1 ^ 
h n = 4 11 



/ 



V 



sin (y^f 1 ) 



sm 



23, 5.4.15.1] we find that Lo(ip } dj) is given by the same expression as in Section 3.1. 



In this way, the sought resolvent kernel is given by a formula analogous to ( |3.2|) with 

X(a,a, <p] z) = A(a, a, 0; z)' 1 , (4.6) 
where Aj r = — Go(dj, a r ) for j ^ r , while the diagonal elements are given by 



Ajj = aty - lim (G (x, dj) + — In \x - aj 



The limit is easily evaluated using (|4.5| ); we arrive at the expressions 

sinh(27TA/^ 2 ^ 2 — z) sin 2 (vnbj) 



n=l 



^cosh(27Tv / ' /2 ^ 2 — z) — cos0 d\J u 2 n 2 — z 2ixn J 

(4.7) 

A jr = -G (d j} d r ) , j ^ r ; 

notice that the matrix A(z) does not change if the flux <p is modified by an integer 
multiple of 2ir . 



4.3 Spectral properties 

Since the point interactions represent a finite-rank perturbation to the resolvent, 
a ess (H(a, a, (p)) = a ess (Ho((j))) = , i.e., the spectrum of if (a, a, 0) is again purely 
discrete. 

Moreover, it is disjoint with 8 := { e mn : m E n — 1, 2, . . . } provided there 
is at least one point perturbation present. If some e mn would be an eigenvalue 
of H(a, a, 0) , the corresponding eigenfunction had to satisfy the corresponding 
Helmholz equation everywhere outside the set {a±, . . . , a^} being therefore a finite 
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linear combination tp = J2m,n c mnVm ® Xn ■ Such a function, however, is bounded in 
Q and cannot have a logarithmic singularity, i.e., L o ('0,a 3 ) = for any j . 

A real number z ^ £ therefore belongs to the resolvent set of H(a, a, 0) , so we 
are allowed to use the argument of Section 3.3. The eigenvalues are given by the 
condition 

det A(a,a, <p;z) = 0, (4.8) 

under which the linear system 

N 

A{z) jr d r = 

r=l 

has a solution (or solutions); the corresponding (non-normalized) eigenfunctions are 
then 

N 

<Pffl = J2 djG (x,aj\<l>,z) . (4.9) 

r=l 

5 Point perturbations on an infinite cylinder 

If the cylinder heigth d of the previous section goes to infinity, the configuration 
space becomes again noncompact; one can map it on the strip 

Q = [0,2tt) x M. 

The results of Section |^ modify easily to the present situation. 

5.1 The resolvent and scattering amplitudes 

The free Hamiltonian is defined again by the cyclic boundary conditions ( |4.2| ), so 
the "transverse eigenfunctions" are r) m : r) m (x) = (2'K)~ l l 2 e^ m+A > x corresponding 
to the eigenvalues e m := (m + A) 2 , m G % . The free resolvent kernel is therefore 

j e i{m+A)(xi—X2) e ikm(z)\yi— j/a| 

where 



k m (z) := yjz-{m + A) 2 . (5.2) 

Introducing now iV point perturbations as in Section |3|, we find that the full resol- 
vent kernel is given by the relations ( |3.2|) and Q3.4|) with 

. i 



1 / 1 J_ 

k o(z) \k m (z) \m\ 



(5.3) 

A jk = -G Q {a,j,a k ) , j ^ k ; 
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Figure 9: A strip potential due to an external field. 



It is convenient for computational purposes to rewrite the first of these relations in 
the form 



i 

Ait 



k (z) 



E 



k m (z) \m+A\ 



where otj are renormalized coupling constants, 



(5.4) 



«,■ = OLi 



+ S E Ira 



A 



\m+A\ 



we suppose, of course, that A £ 2fi . 

From here one can find the discrete spectrum of H(a, a, <fi) by virtue of the eqs. 
and (|4.9|) . Similarly, the reflection and transmission amplitudes are 



N 



"Tn.m. \ Z ) 



tnm (Z) 



Air 



E (A^r 1 ), 



ik~ 



j,k=l 



i(m+A)bj gi(n+A)fej. 
km. (z) 



(5.5) 



N 



5nm + T= E ( A (^) ")jk- 
41T j,k=l 



the elementary block S nm of the S-matrix is given again by ( p.24| ), where the tilded 
quantities are obtained by mirror transformation. The scattering is, of course, a 
one-dimensional problem with an additional degree of freedom corresponding to the 
transverse channels — it imposes no restriction on the size of the cylinder contrary 
to the claim made in 



6 An infinite strip in external field 

Let us return now to the infinite strip and assume that, in addition to point per- 
turbations, the particle interacts with a potential. Such a situation can occur, for 
instance, if the strip is cut not from a plane but from a curved surface (see Fig. 9) 
and exposed to a homogeneous electric field. We shall suppose that the potential 
depends only on the longitudinal variable; otherwise the problem would become 
much more complicated. 

6.1 The unperturbed Hamiltonian 

In the described situation the Hamiltonian free of point perturbations is therefore 

H v := -d 2 x -d 2 y + V{x) (6.1) 



2s 



with the usual Dirichlet condition (|2.1| ) at the strip boundary. Our assumptions 
about the potential are for the moment rather weak: we require only that V G 
Lj oc (]R) and the operator 

h v ■= -d 2 x + V{x) 

with the natural domain is e.s.a. Then Hy understood as the closure of hy(&(— d 2 ) 
is self-adjoint and its spectrum is 

a(H v ) = {X + n 2 : A G a(h v ) , n = 1, 2, . . . } . (6.2) 

In particular, the resolvent set p(H v ) D (— oo, V- + 1) , where V_ := mf xe]R V(x) . 
Of course, the character of the spectrum depends strongly on that of a(Hy) . If, for 
example, the potential vanishes at infinity, lim[ x [_+ D0 V(x) = , then hy is below 
bounded with a ess (hy) = [0, oo) by |25|, Sees. X. 2 and XIII. 4], so 

a ess (H v ) = [1, oo) 

as in the absence of the potential. In a similar way, if the potential grows at large 
distances (at least in the mean) then the spectrum of Hy is pure point etc. 

By assumption, hy is limit point at both ±oo , so it follows from the Weyl 
alternative that to any z = k 2 G € \ IR there are functions u = u(-,k) and 
v = v (•, k) , unique up to a multiplicative constant, which solve the equation /" + 
{k 2 — v)f = being square integrable at =Foo , respectively. They are pointwise 
continuous with respect to k ; we use the same symbols for their extensions to the 
positive real axis in the /c-plane. We denote k n (z) := y/z — n 2 as before and set 

u n := u(-,k n (z)), v n := v(-,k n (z)). 

The free resolvent kernel is then 

G v {xi,x 2 ;z) = -- V Un &<) Vn ( x >) sin(nyi) sin(y 2 ) , (6.3) 

7T ^ W{U n ,V n ) 

where the smaller and larger of xi, x% , respectively. 

6.2 Point perturbations 

The Hamiltonian H(a, a, V) with point perturbations is defined again through the 
boundary conditions ( [3.1|) . To find its resolvent, we need the logarithmic singularity 
of the kernel ( |S.3| ). If the potential is smooth the latter is known to be the same 
as in the absence of the potential |p9[ . In fact the potential need not be regular; 
nevertheless an additional restriction is required because the property under con- 
sideration is local, and therefore could be affected by a local singularity. We shall 
assume that the potential is locally bounded, i.e., bounded in any compact interval. 

To find the singularity, we notice that what matters here is the behavior of the 
series at large n ; however, these terms correspond to solutions which are deep in the 
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classically forbidden region (i.e., with the imaginary part of the channel momentum 
large positive) and therefore of the WKB-form, 

Un(x) ^(sj-^er^'Erny, v n (x) = ^(x)' 1 / 2 e ~ JT Err^/c) , 

where n n (x) := \Jn 2 + V^x) — z and Err Q (/t) stands as a shorthand for l + , 
so 

u n (x)v n {x) 1 1 

"7777 T = ~o 7-T Err l («n) = - — Err l n 

due to the local boundedness of the potential. Hence we obtain 

G v (x,dj;z) 1 
hm — ^— - = , 

|x-a,-|-t-0 m(x — aj| 27T 

and the argument of Section 3.1 can be easily modified giving 

N 

(H(a, a, y)-z) _1 (fi,f 2 ) = G v (x 1 ,x 2 ;z)+ ^ (A(z) _1 ) jfe GV(£i, a,-; z)G v (a k , x 2 ; z) , 

i,k=i 

(6.4) 

where A (a, a, V; 2) is the following N x N matrix: 

A jr = -G y (a i ,a r ;z) , j ^ r , 



(6.5) 



It becomes singular if the condition 

det A(a, a, V; z) = (6.6) 

is valid; the corresponding z are eigenvalues of H(a, a, V) . If z £ p(H v ) , the 
argument of Section 3.2 may be used: the non-normalized eigenfunctions are of the 
form 

N 

= ^djGvix^^z) , (6.7) 
where d £ G N is the respective solution to J2k=i ^( z )jkdk = . 
6.3 Scattering 

One has to assume that the underlying scattering problem for the potential V 
makes sense, the wave operators exist and are complete; since V is locally bounded 
by assumption, this is ensured, for instance, if 



\V(x)\ < C(l + |x|)- 



-l-e 
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for some positive C, e |25, Sec. XI. 4]. Repeating the limiting-absorption argument 
of Sees. 2. 6 and 3.4, we get the generalized eigenfunction 



N 

^(f) = e**W-x»(y) + E (A(z)- 1 ) jk G v (x,a j ;z)e lk ^ X n(b k ). 

j,k=i 

We substitute for Gy from ( |6.3| ) and consider the asymptotic behavior. To this 
end, we need the asymptotics of the elementary solutions which is given by 



. g ik n (z)x ^ 

U n (x) « s| T ^ e -ik n (z)x _^_-p e ik n (z)x x _^ 



and an analogous formula with the interchange of signs for v n (x) ; here r n : = 
r(k n (z)) and p n := p{k n {z)) are the standard (left-to-right) transmission and re- 
flection amplitudes for the potential V at the momentum k n (z) . This yields 

r (z) = - T(Mz)- 1 U ^^ i) ^ nbk K (a-)e ik ^ 
nm{ ) ^ 3 h hk k m (z)r(k m (z)) mK l) 

(6.9) 

± / x r i v^/a/\-i\ sin(mfo 7 ) sin(nfo fc ) . . ; fr 
vr ^ x k m {z)T{k m (z)) 



the unitarity condition ( p. 20 ) remains the same. 



7 Periodic perturbations 

An extension of the results discussed above to an infinite number of perturbations is 
straightforward for the algebraic part; the matrices A become infinite-dimensional, 
i.e., operators on the Hilbert space £ 2 . However, since they are generally unbounded, 
the corresponding analysis is much harder. The problem should be treated along the 
same lines as in |], Sec. III. 1.1] but we are not going to do that; instead we concentrate 
on an important particular case of a strip in which the point perturbations are 
periodically arranged. 

7.1 The Floquet— Bloch analysis 

Suppose that the set of point perturbations {a, a} = {[%, Uj] : j = 1,2,...} in fl 
is countably infinite and has a periodic pattern with a period £ > . As usually we 
suppose that the width of the strip is d = ir ; the general case can be obtained by 
a simultaneous rescaling of the coupling constants. 

The existence of the corresponding operator H(a, a) can be established as 
in Sec. III. 1.1]. Following the standard procedure named Floquet by mathe- 
maticians and Bloch by physicists, we can find a unitary operator U : L 2 (Q) — ► 
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L 2 (B, (£/2n)d9; L 2 (ft)) , where 

Cl ■= [0,£] x [0,tt] , B :-- 



the x-projections of these rectangles are the Wigner-Seitz cell of the underlying 
one-dimensional lattice, and the corresponding Brillouin zone, respectively — cf. 
|], Sec. III. 1.3]. By means of U , the operator H(a,a) is unitarily equivalent to 



UHia.ajU- 1 = — / H(a,a;9)d9, (7.2) 

27T J\9£\<n 

where H(a, a; 9) is the point-interaction Hamiltonian on L 2 {Vt) , i.e., the Laplacian 



satisfying the boundary conditions ( |2.2j ) for a: G [0 



i>(£-y) = e m ^(0+,y), ^(f-v) = e m ^(0+,y) (7.3) 

for y G [0, 7r] and (|3.1|) ; with an abuse of notation we shall use the symbol (a, a) 
for the subset of the point interactions in Cl and suppose that their number is iV . 

The operators H(a, a; 9) are easy to investigate because up to a change in no- 
tation and a simple rescaling they coincide with the magnetic cylinder Hamiltonians 
of Section 4. In particular, the "free" eigenvalues 



/ 27rm \ 2 

e m n(0) : = (— + 9) + n 2 , meZ,n=l,2,..., (7.4) 

correspond to the eigenfunctions rj e m ®Xm where \n are the usual transverse solu- 
tions and 

rj° m (x) := J-eiptm+oWt, m G Z . 
V £ 

Moreover, the free resolvent kernel is 

sinh((£— \x\ — X2\)\/ri 2 — z) + e 2ir,de sinh( \x\ — x 2 \ Vn 2 — z) 



Gq{x x ,x 2 ]9]z) = 



,. :] cosh(^V^ 2 — z) — cos( 

(7.5) 

sin(m/i) sin(nj/2) 



x 



iiy/n 2 — z 



where r\ := sgn {x\— x<i) , and the full kernel is expressed by a formula analogous to 
(|3) with 

X(a,a, 9;z) = A(a,a,9; z)^ 1 , (7.6) 
where A Jr = — G Q (dj, a r ) for j ^ r , while the diagonal elements are given by 

1 ~ / sinh( l\[n 2 ^~z) sin 2 n6 ? - l\ 

A - • = a , V / '- 3 . (7.7) 

" 3 vr ^ V cosh (iVrf^ -cos 9£ Vri^ 2n y ' 
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We may also write the last formula as Ajj(a, 3, 0; z) = ctj — £(3j, 0; z) , where the 
function £ is more explicitly given by 



£(3*,0;z) = - 



1 ^ / sin(£y/ z — n 2 ) sin 2 nbj 1 



» " V cos(£y/ z — n 2 ) - cos0£ v 'z-ri 2 2n / 



n=l 



(7A 



1 ~ / sinh(£ v / r^I) sin 2 nbj 1 N 
TT n= fd +1 Vcosh(V« 2 --2) - cos0£ V^ 2 --2 2n / 



n=n[.z] + l 

we have denoted here n[z] := min{0, • It is defined everywhere except at 

£(3,0) := {e mn (0)G£(0) : sinn&^O}, (7.9) 

where £(0) is the eigenvalue set Q7.4Q ; it is also easy to see that the function £ is 
monotonously increasing between any pair of neighboring singularities. 

To find the spectrum of the original operator, one has to solve first the spectral 
problem for H(a, 3; 0) in analogy with the magnetic cylinder case: the eigenvalues 
are given by 

det A(a,a,6\z) = 0, (7.10) 

and the corresponding eigenfunctions are given by a relation similar to (|4.9|) . Next 
one has to overlay the obtained eigenvalue sets for running through the x- 
projection of the Brillouin zone to get the spectrum of H(a, 3) . The main question 
is, of course, whether we cover a single interval in this way as in the free case, or 
whether the interaction will open some gaps. 



7.2 A single array of perturbations 

Since it is not easy to draw general conclusions about the form of the spectrum, let 
us restrict to the simplest case where iV = 1, i.e., the strip Q contains a single 
periodic array of point perturbation of a strength a and spacing i . The condition 
Q7.10D then simplifies to 

£{a,6;z) = a. (7.11) 

We have seen that the function on the left hand side is monotonously increasing 
between its singularities, i.e., the points of £(3,0) . This means that for fixed a, 
there is a sequence {e r (a, a, 0)},^ o arranged in the ascending order; each of them 
depends, in fact, only on the y-component b of the vector 3. The lowest one 
satisfies 

e o (a,a,0) < 1 + 2 , 

and between each two neighboring points of £(3,0) there is just one of the other 
eigenvalues. It is also clear that any of e r (a, a, 0) is continuous with respect to the 
parameters and e r (-,3, 0) is increasing for fixed b and 0. 



33 



On the other hand, the eigenvalues need not be monotonous with respect to 9 . 
The implicit-function theorem tells us, however, that 

de r (a,a,9) d£(a,9;z) ( <9£(a, 6>; z)\ 1 



d9 d9 \ dz 



(£r,0) 



whenever the denominator is nonzero. Leaving aside the thresholds, z = n 2 , and 
the points of £ (9) , a straightforward differentiation shows that £(o, •;•) is ana- 
lytic in both variables. Since the function £ is not identically zero, the derivative 
de r (a,a,9)/d9 may be zero at some points but never in an interval. This means 
that the spectrum of H(a,a) is absolutely continuous — cf. p5| , Sec. XIII. 16]. 



Let us turn now to the question about the number of gaps. Below z = 1 the 
spectrum may be estimated by means of extrema of the function £ : we have 

\ 1 \ 



£+(a, z) := max £(a, 9: z) = — . coth -\/n 2 — z , . 

and a similar formula for the minimum, £_ (a, z) , with coth replaced by tanh . 
Both functions are increasing and have the same logarithmic asymptotics as £ (a, z) 
of Section 2.5 as z — ► — oo . On the other hand, £+(a, •) diverges as z — > 1— , while 
£_(a, •) has a finite limit. This shows that a gap exists provided 



£ 1 
a < £-(a, 1— ) = — sin 2 6— — 

sin 2 nfe 



(7.12) 




tanh —\/n 2 — 1 — — 

v/^T V 2 / 2 ^/ 



The condition is satisfied for a strong enough coupling, or alternatively, for any fixed 
a and the spacing I large enough. 

This conclusion calls to mind the spectrum of a straight polymer in a plane 
described in Sec. III. 4] which has also one gap if the coupling is stronger than 
some critical value. However, for a "coated polymer", i.e., our strip with an array 
of perturbations a much stronger result is valid; we are going to show that under a 
suitable choice of parameters it can have any finite number of gaps. 

To this end, we consider z G (2, 3) and £ 1 , and rewrite the right hand side 



of the relation ( |7.8|) as 

f (a, 9; z) = f (a, 0; z) + ij(a, 9; z) 

where 

fo(a, 0; z) 



sin(£vz—T) sin b 



cos(£\/z— 1) —co§9£ tx\J z — \ 

and rj(a,9;z) is the rest. In the same way as above, one can show that r)(a,9; ■) 
is monotonously increasing with a bounded derivative everywhere below the second 
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threshold, in particular, in the chosen interval of energies. Moreover, 



^ ■ ■ \ I sin 2 nb . (I r^, \ 1 

rj{a, 6;z) < r] + {z) := - — + - > — coth 




n=2 



and the minimum, rj-(z) , is obtained when coth is replaced by tanh . These 
estimates become close if £ is large; using the inequality cothw — tanhw < 5e~ 2u 
for 2u > 1 , we find 



rj + (z)-V-(z) <~t < ' t e 

n n =2 Vn 2 -z n n=2 



-£(n^l) 



for z G (2, 3) , so 



V+( z ) ~V-( Z ) < 



ti 1 — e 1 

and the allowed corridor shrinks exponentially with increasing I . On the other 
hand, the function 

, . sin u 

ge ^ u > := To 

cos u — COS Ul 

is increasing between any two zeros of its denominator. In the intervals, where it is 
positive, it is estimated by the appropriate branch of tan (| + Trrnj from below; if 
it is negative, we have a similar estimate from above with tan replaced by — cot . 
Hence independently of 6 we have either 

. ^ „ sin 2 b . , 

4o(a, 0; z) > — -j== tan - 

or 

j_ , ^ „ . sin 2 b , 

£o{a,0;z) < -j== cot - 

TTy/z— 1 

where {•} denotes the fractional part. Putting the estimates together, we see that 
the oscillating part dominates, so for sufficiently large \a\ there are gaps having 

1+ (Jy^ as one endpoint provided the latter lies in the chosen interval. In addition, 
tan u + cot u > 2 , which means that the gap between the lower and the upper bound 
to £(a,9;z) never closes for z e (2,3) if 

— 7 < V2 sin 2 b . 




Together we infer that for any a G M, the Hamiltonian under consideration can 
have an arbitrary finite number of gaps in its spectrum provided i is chosen large 
enough. 

Recall that the number of gaps is a quantity strongly dependent on the dimen- 



sion of a periodic system. For one-dimensional systems it is generically infinite [25 



Sec. XIII. 16], |TJ, Sec. III. 2], while for higher dimensions the Bethe-Sommerfeld con- 
jecture states this number is finite. This is known to be true for regular periodic 
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potentials — see, e.g., |27], |28fl . If the potential is replaced by a point interac- 
tion, the one-dimensional situation suggests that the gap closing at high energies is 
slower. The Bethe-Sommerfeld conjecture is known to be valid for two- and three- 
dimensional lattices if each lattice cell contains a single point perturbation — cf. [|l], 
Secs.III.l, III. 4]. On the other hand, an example of a "multidimensional" periodic 
graph system has been presented recently |13[ [14| which exhibits an infinite num- 
ber of gaps — or none at all. In this connection it is interesting to ask whether a 
periodically perturbed strip considered here can exhibit an infinite number of gaps. 
However, we are not going to attack this problem in the present paper. 



8 Random point perturbations 

Again, we have no intention to develop a general theory which would modify the 
results of Sec. III. 5 in Jl]] to the present situation. Our goal is more modest; we 
want to investigate numerically how the scattering changes in the above studied 
models if the set of point perturbations is chosen at random at an interval of a 
length L . We will thus investigate a family of random operators = H(a, a w , (p) , 
where (ai ja; , a 2 ,u,, . . . , a^^) and (bi jU1 , &2,u>, • • • , b^^) are independent random num- 
bers, a,j >u G (0, L) , bj )W G (0, 2tt) for all j = 1, . . . , N . 

It is widely accepted that an irregular scattering on such a family of impurities 
leads to fluctuations of the corresponding S-matrix, and consequently, of the mea- 
sured sample conductance. These mesoscopic fluctuations are generic in the sense 
that they do not depend on detailed properties of the sample under investigation. 
Our aim here is to find the conductance fluctuations for the present model. 

The mentioned universality is usually manifested through relations to spectral 
properties of random-matrix ensembles. More specifically, the S-matrix eigenvalues 
of an irregular quantum scattering system are usually supposed to conform with 
those of a circular ensemble of random matrices. This correspondence has theoreti- 
cally been demonstrated by Blumel and Smilansky for the first time using semiclas- 
sical arguments [0. According to it, for systems with time-reversal symmetry S 
corresponds to the circular orthogonal ensemble (COE), while for systems without 
this symmetry its statistical properties refer to the circular unitary ensemble (CUE). 
However, as we have said, the use of the COE/CUE ensembles for description of the 
mesoscopic fluctuations was justified with the help of a semiclassical approximation 
[]7|, |5|, which implicitly requires a large number of open channels. It is therefore in- 
teresting to investigate the situation, when only a few transverse modes contribute 
to the scattering, so the mentioned argument may not be applied. 

The quantity to study is the sample conductance which is obtained (in the two- 
probe setting) from the transmission amplitudes ( |3.19| ) by the Landauer formula 
( |3.21|) . A standard approach to conductance fluctuations consists of regarding the S- 
matrix as a random matrix of a dimension given by the number of open channels and 
averaging the corresponding conductance over the appropriate random ensemble; 
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Figure 10: var(G) for random perturbations vs. momentum 



Figure 11: Conductance for random perturbations vs. momentum 
this leads to [HI: 



M M 



for COE 



2 4M+2 

(G) = { (8.1) 

for CUE, 



M 
2 



where (G) is the mean value and M is the number of open channels in the cor- 
responding quantum wire. For M — > oo we get (G)cue — (G)coe = j, i-e-, the 
averaged conductance is enhanced by 1/4 , if the time reversal symmetry is broken. 
For the variance var(G) the averaging over the circular ensembles leads to 



M(M+1) 2 
(2M+l) 2 (2M+3) 

var(G) = { (8.2) 



2 



M 

4(4Af2-l) 



for COE 
for CUE. 



This relation shows that, in the semiclassical case, the conductance variance is prac- 
tically independent of M being equal to 1/8 and 1/16 for COE and CUE, respec- 
tively. 

The model we are considering here offers an opportunity for testing these re- 
lations in a situation where the full solution to the quantum scattering problem is 
known and the semiclassical approximation is not required. To be specific, we will 
consider the dependence of the conductance on energy of the incomimg particle in 
the case when the number of impurities and their distribution inside the waveguide 
is fixed; we choose 25 impurities distributed at random inside a waveguide segment 
of length 25. 

One expects that for a fixed number of impurities the universality regime de- 
scribed above can be valid only within a certain energy interval. For energies which 
are too low the localization effects will dominate the scattering scenario and the 
incoming wave will be nearly completely reflected. On the other hand for energies 
above a certain threshold direct processes will dominate and the major part of the 
incident wave function will be transmitted through the sample; the region interme- 
diate between these two regimes is expected to be the universality "window" . 

Our numerical results support the above conjectured picture. In Figure 10 we 
plot var(G) as a function of the incident wave vector k in comparison with the 
random-matrix prediction (|8.2j ). The universality window is localized approximately 
between k = 3.5 and k = 5.5. For smaller k the localization effects suppress the 
variance whereas for k above the window the variance depends significantly on k 
with maxima localized at integer values of k . 



Figure 12: Level-spacing statistics for random-impurity scattering 
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The conductance G of the same randomly perturbed waveguide is plotted on 
the Figure 11. The prediction of ( |3.6| ) is again shown for the sake of comparison 
(the dashed line). It is clear that G coincides with the random-matrix prediction 
only inside the universality window; for smaller k it is suppressed by localization 
effects while for large k it is enhanced by direct transmission. 

Finally, the level-spacing statistics of the corresponding S-matrix is shown on 
Figure 12. The first picture corresponds to k inside the universality window, while 
the other one is computed for a larger momentum value. The prediction of the 
random-matrix theory (Wigner distribution) and the Poisson distribution are plot- 
ted for the sake of comparison; a shift towards the latter for k above the window is 
clearly visible. 
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Figure captions 

Figure 1 The dependence of £(a, z), where a = (0, b). The dotted, dash-dotted 
and dashed curves represent b = tt/2, b = tt/6 and b = 7r/12, respectively. 
The full line is the leading term of asymptotics, £, as (a, z) = — In |) + 

Figure 2 Eigenvalue plots for different values of b, a is the coupling constant 
- cf. (2.13). The dashed, full and dash-dotted curves represent b = 7r/2, 
b = 7r/6 and b = n/12, respectively. 

Figure 3 Unnormalized eigenfunctions for different values of b, a = 1. From 
top to bottom, they belong to b = tt/2, b = n/6 and b = tc/ 12, respectively. 

Figure 4 Pole trajectory in the z-plane. The perturbation is located in the 
middle of strip, i.e. b = tt/2. The resonance is on the third sheet of Riemann 
surface. The trajectory begins at the treshold z — 9, which corresponds to 
a — > +oo and runs through the lower part of the complex plane as a decreases 
to — oo. In the inset is the dependence of the residue p := £'(a,z) depending 
on the coupling constant a . 

Figure 5 Eigenfunctions for three point perturbations. The ground state, E = 
-8.7051 (bottom) and the second excited state, E 2 = -7.8199 (top). The 
perturbation are placed at (0, 7r/2), (1, 2ir/3) and (2, 7r/4); they have the same 
strength a = —0.15. 
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Figure 6 An embedded eigenvalue due to symmetry. The ground and first 
excited energies of a strip with two perturbations placed symmetrically with 
respect to the axis of strip at (0, 7r/3) and (0,27r/3). As a increases, the first 
excited state reaches the continuum and tends to 4 as a — > oo. 

Figure 7 The cascading effect for a cluster of eigenvalues. Four perturbations 
placed at (— l,27r/3), (0, 7r/2), (l,7r/3) and (2, 27r/5); the first three have the 
same strength oti = a 2 = a 3 = —0.2. If we change a 4 from +oo (very 
weak perturbation) to — oo (strong perturbation), there is a critical value at 
which the third excited states appears below the continuum. It becomes more 
bound as decreases, joins the triplet E , E 1 , E 2 and the ground state energy 
E parts the cluster and drops exponentially. The order of levels is always 
preserved. 

Figure 8 The conductance plots for a pair of point perturbations. The plots 
show the dependence of conductance G on energy z. The perturbations are at 
(0, 7r/3), (0.25, 27r/3), having the same strength a = 0.2 (the upper plot) and 
at (0,6tt/11), (10,37r/ll), ot x = -0.25, a 2 = -0.2 (the lower plot). 

Figure 9 A strip potential due to an external field. 

Figure 10 var(G) for random perturbations vs. momentum. 

Figure 11 Conductance for random perturbations vs. momentum. 

Figure 12 Level-spacing statistics for random impurity-scattering. 



41 



